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Abstract 



intrabeam scattering growth rates be 



Beginning with the general Bjorken-Mtingwa solution, ^ 
we derive a simplified model of intrabeam scattering (IBS), 
one valid for high energy beams in normal storage rings; 
our result is similar, though more accurate than a model due 
to Raubenheimer. In addition, we show that a modified ver- 
sion of Piwinski's IBS formulation (where r/ 2 y / (3 x ,y has 
been replaced by Ti x ,y) at high energies asymptotically ap- 
proaches the same result. 
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with dp the relative energy spread, e x the horizontal emit- 
tance, and e y the vertical emittance. The growth rates ac- 
cording to Bjorken-Mtingwa (including a y2 correction 
factor Jl3|], and including vertical dispersion) are 
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1 INTRODUCTION 

Intrabeam scattering (IBS), an effect that tends to in- TrL^'Tr 
crease the beam emittance, is important in hadroniclQ and 
heavy ion[||] circular machines, as well as in low emittance 
electron storage rings [|3}. In the former type of machines 
it results in emittances that continually increase with time; 
in the latter type, in steady-state emittances that are larger 
than those given by quantum excitation/synchrotron radia- 
tion alone. 

The theory of intrabeam scattering for accelerators was 
first developed by Piwinski[[||, a result that was extended 
by Martini[|], to give a formulation that we call here the 
standard Piwinski (P) method^; this was followed by 
the equally detailed Bjorken and Mtingwa (B-M) result[^|. 
Both approaches solve the local, two-particle Coulomb 
scattering problem for (six-dimensional) Gaussian, uncou- 
pled beams, but the two results appear to be different; of the 
two, the B-M result is thought to be the more general^. 

For both the P and the B-M methods solving for the 
IBS growth rates is time consuming, involving, at each 
time (or iteration) step, a numerical integration at every 
lattice element. Therefore, simpler, more approximate 
formulations of IBS have been developed over the years: 
there are approximate solutions of Parzen [[)]], Le DuffjlO||, 
Raubenheimer JIT) , and Wei[|l2|. In the present report we 
derive — starting with the general B-M formalism — another 
approximation, one valid for high energy beams and more 
accurate than Raubenheimer' s approximation. We, in addi- 
tion, demonstrate that under these same conditions a mod- 
ified version of Piwinski's IBS formulation asymptotically 
becomes equal to this result. 

2 HIGH ENERGY APPROXIMATION 

2.1 The General B-M Solution^ 

Let us consider bunched beams that are uncoupled, and 
include vertical dispersion due to e.g. orbit errors. Let the fj^ _|_ y ^ 
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where i represents p, x, or y; 
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with ro the classical particle radius, c the speed of light, N 
the bunch population, (3 the velocity over c, 7 the Lorentz 
energy factor, and er s the bunch length; (log) represents 
the Coulomb log factor, () means that the enclosed quan- 
tities, combinations of beam parameters and lattice prop- 
erties, are averaged around the entire ring; det and Tr sig- 
nify, respectively, the determinant and the trace of a matrix, 
and I is the unit matrix. Auxiliary matrices are defined as 
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The dispersion invariant is H = [rj 2 + (/??/ — ^(3'r]) 2 }/ f3, 
and (j> = f]' — ^P'ij/P, where (3 and r\ are the beta and 
dispersion lattice functions. 

The Bjorken-Mtingwa Solution at High Energies 

Let us first consider 1/T V as given by Eq. ^| Note that if 
we change the integration variable to A' = Xa 2 H / 7 2 then 

a 2 + X' ~a( x 
-aQ x 1 + A' -bC, y 
-bC,y b 2 + A' 



H 



(8) 



2 





x,y 



V x iV 

(10) 

Note that, other than a multiplicative factor, there are only 
4 parameters in this matrix: a, b, £ x , ( x . Note that, since 
[3cj) 2 < H, the parameters ( < 1; and that if TL k, r) 2 //3 
then C is small. We give, in Table 1, average values of a, b, 
( x , in selected electron rings. 

Table 1: Average values of a, b, ( x , in selected electron 
rings. The zero current emittance ratio ~ 0.5% in all cases. 



Machine 


E[GeV] 


iV[10 lu ] 


(a) 


(b) 


<Q 


KEK's ATF 


1.4 


.9 


.01 


.10 


.15 


NLC 


2.0 


.75 


.01 


.20 


.40 


ALS 


1.0 


5. 


.015 


.25 


.15 



Let us limit consideration to high energies, specifically 
let us assume a,b <C 1 (if the beam is cooler longitudinally 
than transversely, then this is satisfied). We note that all 
3 rings in Table 1, on average, satisfy this condition rea- 
sonably well. Assuming this condition, the 2nd term in 
the braces of Eq. ^| is small compared to the first term, and 
we drop it. Our second assumption is to drop off-diagonal 
terms (let ( = 0), and then all matrices will be diagonal. 

Simplifying the remaining integral by applying the high 
energy assumption we finally obtain 
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A plot of g(a) over the interval [0 < a < 1] is given in 
Fig. U to obtain the results for a > 1, note that g(a) = 
g(l/a). A fit to g, 
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is given by the dashes in Fig. [|. The fit has a maximum 
error of 1.5% over [0.02 < a < 1]. 

Similarly, beginning with the 2nd and 3rd of Eqs. ^, we 
obtain 



1 _ <T 2 p (H x ,y) 1 
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Our approximate IBS solution is Eqs. |ll 14 Note that 
Parzen's high energy formula is a similar, though more ap- 
proximate, result to that given here|]]; and Raubenheimer's 
approximation is Eq. [n], with g(a/b)aH )& p replaced by i, 
and Eqs. 14 exactly as given here [|TT]] . 

Note that the beam properties in Eqs. [TT]JT4j , need to be 
the self-consistent values. Thus, for example, to find the 




Figure 1: The auxiliary function g(a) (solid curve) and 
the approximation, g = a ( 021 -° 0441nQ ) (dashes). 

steady-state growth rates in electron machines, iteration 
will be required^]. Note also that these equations assume 
that the zero-current vertical emittance is due mainly to ver- 
tical dispersion caused by orbit errors; if it is due mainly to 
(weak) x-y coupling we let Tiy — 0, drop the 1 /T y equa- 
tion, and let e y — ne x , with n the coupling factor^. 

What sort of error does our model produce? Consider a 
position in the ring where ( y — 0. In Fig. || we plot the 
ratio of the local growth rate T^ 1 as given by our model to 
that given by Eq. ^ as function of for example combi- 
nations of a and b. We see that for ( x < Vbe*- 1 -^ (which 
is typically true in storage rings) the dependance on ( x is 
weak and can be ignored. In this region we see that the 
model approaches B-M from above as a,b — > 0. Finally, 
adding small C, y ^ will reduce slightly the ratio of Fig. 0. 
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Figure 2: The ratio of local growth rates in p as function 
of Cx, for b = 0.1 (blue) and b = 0.2 (red) [( y = 0]. 

3 COMPARISON TO PIWINSKI 

3. 1 The Standard Piwinski Solution^ 
The standard Piwinski solution is 
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the function / is given by: 

f(a,b,q) = 8tt / riu 1 ^ x 
Jo r V 
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The parameter d functions as a maximum impact parame- 
ter, and is normally taken as the vertical beam size. 

3.2 Comparison of Modified Piwinski to the B- 
M Solution at High Energies 

We note that Piwinski's result depends on rj 2 / (3, and not 
on H and <f>, as the B-M result does. This may suffice for 
rings with (7i) w (rj 1 / j3). For a general comparison, how- 
ever, let us consider a formulation that we call the modified 
Piwinski solution. It is the standard version of Piwinski, 
but with rj 2 /P replaced by H (i.e. a, b, ah, become a, b, 
<jh, respectively). 

Let us consider high energy beams, i.e. let a,b <C 1: 
First, notice that in the integral of the auxiliary function / 
(Eq. the -0.577 can be replaced by 0; the -3u 2 in 
the numerator can be set to 0; P (Q) can be replaced by 
y/a 2 + u 2 (vF + t?). The first term in the braces can be 
approximated by a constant and then be pulled out of the 
integral; it becomes the effective Coulomb log factor. Note 
that for the proper choice of the Piwinski parameter d, the 
effective Coulomb log can be made the same as the B-M 
parameter (log). For flat beams (a <C b), the Coulomb log 
of Piwinski becomes (log) = In [da 2 H / (4r a 2 )]. 

We finally obtain, for the first of Eqs. 
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We see that the the approximate equation for 1/T p for high 
energy beams according to modified Piwinski is the same 
as that for B-M, except that h(a, b) replaces g(a/b). But 
for a,b small, h(a,b) ~ g(a/b), and the Piwinski result 
approaches the B-M result. For example, for the ATF with 
e v /e x - 0.01, a - 0.01, a/b - 0.1, and h(a, b)/g(a/b) = 
0.97; the agreement is quite good. 



Finally, for the relation between the transverse to longi- 
tudinal growth rates according to modified Piwinski: note 
that for non-zero vertical dispersion the second term in the 
brackets of Eqs. [D] (but with yy 2 y//3 x , y replaced by H. x , y ), 
will tend to dominate over the first term, and the results 
become the same as for the B-M method. 

In summary, we have shown that for high energy beams 
(a,b <C 1), in normal rings (£ not very close to 1): if the 
parameter d in P is chosen to give the same equivalent 
Coulomb log as in B-M, then the modified Piwinski solu- 
tion agrees with the Bjorken-Mtingwa solution. 

4 NUMERICAL COMPARISON® 

We consider a numerical comparison between results of 
the general B-M method, the modified Piwinski method, 
and Eqs. [IT 14 The example is the ATF ring with no cou- 
pling; to generate vertical errors, magnets were randomly 
offset by 15 /im, and the closed orbit was found. For this 
example (Ti y ) = 17 /im, yielding a zero-current emittance 
ratio of 0.7%; the beam current is 3.1 mA. The steady-state 
growth rates according to the 3 methods are given in Ta- 
ble 2. We note that the Piwinski results are 4.5% low, and 
the results of Eqs. [TT[ , p^4j , agree very well with those of B- 
M. Additionally, note that, not only the (averaged) growth 
rates, but even the local growth rates around the ring agree 
well for the three cases. Finally, note that for coupling 
dominated NLC, ALS examples (k — 0.5%, see Table 1) 
the error in the steady-state growth rates (T" 1 ,^ 1 ) ob- 
tained with the model is (12%, 2%), (7%,0%), respectively. 

Table 2: Steady-state IBS growth rates (in [s -1 ]) for an 
ATF example with vertical dispersion due to random errors. 



Method 


l/T P 


1/T X 




Modified Piwinski 


25.9 


24.7 


18.5 


Bjorken-Mtingwa 


27.0 


26.0 


19.4 


Eqs. (TTp4] 


27.4 


26.0 


19.4 
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